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ABSTRACT 

The finding of residual gas in the large central cavity of the HD 142527 disk motivates questions 
on the origin of its non-Keplerian kinematics, and possible connections with planet formation. We 
aim to understand the physical structure that underlies the intra-cavity gaseous flows, guided by new 
molecular-line data in CO(6-5) with unprecedented angular resolutions. Given the warped structure 
inferred from the identification of scattered-light shadows cast on the outer disk, the kinematics are 
consistent, to first order, with axisymmetric accretion onto the inner disk occurring at all azimuth. 

A steady-state accretion profile, fixed at the stellar accretion rate, explains the depth of the cavity 
as traced in CO isotopologues. The abrupt warp and evidence for near free-fall radial flows in HD 
142527 resemble theoretical models for disk tearing, which could be driven by the reported low mass 
companion, whose orbit may be contained in the plane of the inner disk. The companion’s high 
inclination with respect to the massive outer disk could drive Kozai oscillations over long time-scales; 
high-eccentricity periods may perhaps account for the large cavity. While shadowing by the tilted disk 
could imprint an azimuthal modulation in the molecular-line maps, further observations are required 
to ascertain the significance of azimuthal structure in the density field inside the cavity of HD 142527. 

Subject headings: Protoplanetary disks — Planet-disk interactions — Stars: individual: (HD 142527) 


1. INTRODUCTION 

While in its first year of Early Science operations, 
the Atacama Large Millimeter Array (ALMA) already 
provided resolved sub-mm observations of circumstellar 
disks around young stars, where planet formation must 
occur. The pre-main-sequence star HD 142527 (F6llle, 
at a distance of ^145 pc) is particularly interesting as 
its disk is viewed close to face-on, and its protoplanetary 
cavity is amenable to s crutiny thanks to its record size 
(iFukagawa et al.||2QQ6 a|). Indeed, ALMA Cycle 0 obser 
vatio ns of HD 142527 ( [(Jasassus ( 


et al.|2Q13c Perez et al 


2015) revealed gaseous hows inside the dust-depleted cen¬ 
tral cavity, which is a n expected feature of the formation 
of planetary sy stems ( Zhn et a'r]|2Qll Dodson-Robinson 


& Salyk 
with hig 


2011). Here we complete this ALIVIA program 
1 -frequency molecular-line data that provide ad¬ 


ditional insight on the physical structures giving rise to 
the intra-cavity flows in HD 142527. 

An obvious feature of the first ALMA observations of 
HD 142527 is the strong and centrally peaked HCO+(4- 
3) emission, which is consist ent with radial flows in addi- 
tion to Keplerian rotation (Casassus et al.||2013c). In 
this scenario the mass transfer rate inferred frorn the 
HCO+ kinematics matches the stellar accretion rate. It is 
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tempting to associate the HCO+ kinematics with planet- 
induced accretion flows, i.e. gaseous streamers due to 
as yet undetected protoplanets inside the cavity. Such 
accretion kinematics should have counterparts in other 
molecular tracers, and indeed the velocity cent roid map 


in CO(3-2) also hints at non-Kepl erian flows (Casassus 
et al.|2Q13c| ^ iRosenfeld 


d et al.|2Q14 ). However, 


t he ( 


Cycle U beam, in CO(3-2), along with absorption by a dif¬ 
fuse CO screen along the line of sight to HD 142527, pre¬ 
vent adequate sampling of intra-cavity kinematics. Addi¬ 
tionally a warped disk instead of a stellocentric velocity 
compo nent could result in the same CO(3-2) centroid 
maps (Rosenfeld et al. IMit . 

The new ALIVIA observations we present here, acquired 
in the CO(6-5) line, together with a well-determined in¬ 
ner disk orientation thank s to the identificat ion of its 
shadows in scattered-light ( Marino et al.[|2Q15 ), allow us 
to conclude that the intra-cavity hows do correspond to 
stellocentric accretion occurring at all azimuth, and near 
free-fall velocities. Section presents the observational 
results on the intra-cavity gas kinematics. Sec. intro¬ 
duces the parametric models we use to interpret the data, 
and compares the models with the observations. Sec. [4] 
discusses implications of the inferred kinematics in the 
light of the recent theore tical discovery of “ disk tearing” 
in strongly warped disks (Nixon et al.(2013). Sec. [^sum¬ 
marizes our conclusions on the kinematics of stellocentric 
accretion in the cavity of HD 142527. 


2. OBSERVATIONS 

HD 142527 was observed with ALMA in band 9, at 
^700 GHz, during June 2012 and with baselines rang ing 
from 0.04 MA to 0.9 MA. We refer to Appendix | A. 1| for 
an account of data reduction. 
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2.1. Choice of image-plane representations of the 
visibility data 

Image synthesis was performed with a non-parametric 
least-squares fit to the visibility data (see Appendix [B|) ; 
the resulting model images are labelled ‘MEM’. These 
‘deconvolved’ images pick up finer details than the con¬ 
ventional ‘restored’ maps, which are produced by con¬ 
volving the MEM maps with an ellipical Gaussian called 
‘clean’ beam, which is fit to the ‘dirty’ beam (or point- 
spread-function, PSE). With the addition of a dirty map 
of the model residuals, the restored maps allow a stan¬ 
dard representation of the data, with a single and well 
defined angular resolution. However, the moments ex¬ 
tracted on the ‘restored’ maps (Appendixare affected 
by convolution with the clean beam, whicn enhances ve¬ 
locity dispersions by aligning different velocity compo¬ 
nents, with consequences on the intensity maps if ex- 
tractin g mo ments with Gaussian spectral fits (see Ap¬ 
pendix ]b 3|). We thus chose to base our study on the 


moments maps inferred from deconvolved (MEM) dat- 
acubes. A drawback of this choice is that since we are 
extrapolating power in the Eourier domain, the resolu¬ 
tion of the MEM images is variable and depends on the 
image. Here the full-width at half-maximum (EWHM) of 
the MEM PSE varies from 0.062" x 0.051" for a bright 
point source, up to ^50% coarser for fainter extended 
emission (approximated by a collection of spikes). 

2.2. CO(6-5) kinematics 

The main properties of the continuum-subtracted 
GO(6-5) data are summarised by the moment maps 
shown in Eig. The GO(6-5) line intensity displays 
a central decrement so that the minimum intensity is a 
factor of ^ 0.25 ± 0.03 times lower than the average in¬ 
tensity in the brighter ring at a radius of 0.2". This cen¬ 
tral decrement is probably related to photodissociation 
of CO, so that the very central CO(6-5) could be rar¬ 
efied and broad enough to be optically thin. The more 
extended emission in CO(6-5) seems to be ascribed to 
the dust-depleted cavity, with a sharp drop at arcsec, 
corresponding to the radius of the Bright crescent-shaped 
outer ring seen in the continuum. This sharp drop is ab¬ 
sent in CO(3-2), which extends well outside the cavity, so 
that the CO temperature outside the cavity is probably 
not high enough to excite J = 6. 

The disk position angle (PA) is defined as the inter¬ 
section of the plane of a disk with the plane of the sky. 
Eor Keplerian rotation in disks with axial symmetry, the 
velocity centroid maps bear reflection symmetry about 
the disk PA, and the highest Doppler-shifted emission 
towards the blue and red are aligned along the PA. In 
Eig. the velocity centroid is more complex than the 
usuarbutterfly pattern. It is as if the velocity field had 
been ‘twisted’ anti-clockwise, as would happen if the disk 
PA rotated with stellocentric radius r, from —20° East of 
North outside r ^ 0.6", to ^ +60° at r ^0.25", coinci¬ 
dent with the extremes in velocity centroid. Inside 0.1", 
the PA appears to rotate further and reaches 110 ± 10° . 
Yet the scattered light s hadows from the inn er disk show 
that its PA is —8 ± 5° (jMarino et al.||2015). 

The CO(6-5) velocity dispersion reaches very high val¬ 
ues, especially near the star. Velocity dispersion values 
above 4 km s“^ are found in the unresolved central beam; 


the bulk of this dispersion is probably due to flow kine¬ 
matics. However, values of about 2 km s~^ are also found 
away from the central beam, which persist after decon¬ 
volution. Eor instance, in Eig. the 1 a velocity disper¬ 
sion in the restored maps, drest.^ measured at an offset 
P = (+0.2",+0.2"), is 2.01 km s“^, while the 1 a dis¬ 
persion extracted from the deconvolved datacubes, (Jmem, 
still shows a dispersion of 1.78 km s“^. Eurther tests on 
the suitability of the deconvolved images as measures of 
the intrinsic line-width are provided in Appendix [C) 

Beyond 0.1" from the nominal stellar position, the 
line profiles appear to be singly peaked, as illustrated 
in Eig. If due to thermal broadening, the CO velocity 
distribution corresponds to temperatures T = 2Smpa‘^ jk 
in excess of 5000 K, as far out as 0.5" from the star. An¬ 
other mechanism for producing the observed line-widths 
is microscopic turbulence. There are however hints of 
line asymmetry, such as shoulders or lack of low-level 
Gaussian wings, which suggest the presence of multiple 
velocity flows stratified along the line of sight. How¬ 
ever the significance of these spectral features is not high 
enough to motivate detailed interpretation. 

Integrated line profiles for CO(6-5) are shown in Eig.[^ 
The total line flux extracted from the restored datacube, 
using natural weights, and in a circular aperture with 
a radius of 2.5", is 70 Jy km s“^. This value is fairly 
close to the flux densit y of 80 Jy km s“^ reported by 
Casassus et al. (2013b) from APEX+CHAMP+ single- 
dish observations! I'he double-peaked CO(6-5) profiles 
obtained with ALMA are typical of Keplerian disks, and 
are also reminiscent of the CHAMP+ spectrum, given its 
noise level. However, here the blue-shifted peak appears 
to be fainter than the red-shifted peak by ^20%. This 
line asymmetry is opposite to that seen in ^^ CO(2-l) 
(iPerez et al. 2015), but is similar to HCO^( 4-3) dCasas- 
sus et al.||20i3c[ ) and the CO isotopologues ( [Perez et al. 
2015|). It is likely that the lower-J ^^CO is absorbed 
by the diffus e screen at red-shifted velocities (Casassus 
et al. 2013b). The origin of the blue-shifted decrement 


m species not affected by the diffuse molecular screen 
is probably related to the anti-correlation with the con¬ 
tinuum: in continuum-subtracted maps, the line emis¬ 
sion appear s to be quenched under the crescent-shaped 
continuum (Casassus et al. 2013c van der Plas et al. 


2014 
to b 


Perez et al. 


ue-shifted ve 


by Casassus et al 




at locations that correspond 
'deities. The pointing offset invoked 


2013b) to account for the asymmet¬ 


ric TIH7fMI^T"^rdhIeTsprobably unnecessary given the 
match in line flux with the interferometer CO(6-5) maps, 
and given the later finding of the molecular decrements 
at blue-shifted velocities. 


2.3. Complementary sub-mm molecular-line data 

As mentioned in Sec. the HCO+(4-3) and CO(3-2) 
data previously describea in the literature also display 
deviations from Keplerian rotation with a constant disk 
orientation. Although the outer disk follows Keplerian 
rotation with a PA of —20°, inside the cavity the PA join¬ 
ing regions of opposite velocitie^ rotates with stellocen- 
tric distance resulting in ‘twisted’ kinematics (Casassus 


^ relative to systemic 
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Figure 1. Moment maps in CO(6-5) emission, x— and y—axis show angular offset along RA and DEC from the stellar position, in arcsec. 
The moments were extracted on deconvolved MEM models. From left to right we show frequency-integrated line intensity, uncertainty on 
the intensity field, velocity centroid, and 1 a velocity dispersion. 


et al. 2013aj). The high velocit^p] blue and red -shifted 
are offset by 0.2 arcsec, witn a PA of ^60° (Casas- 
sus et al.]|2Q13c ), while the outer disk PA is at ^ —20''. 
Low-velocity filamentary HCO+(4-3), within 0.8 km s“^ 
of systemic, connects the outer disk to the East with 
the central broad-velocity HCO+ (|Casassus et al.|2013c), 
along a PA that coincides with that of the blue-shifted 
signal. A similar but fainter filament may connect with 
the Western outer disk. Interestingly the peak CO(3-2) 
intensity is offset from the star, again towards the same 
direction as the blue-shifted HCO+. The PA joining the 
blue- and red-shifted velocities is thus ^60° and common 
for CO(3-2) and HCO+, and approximately orthogonal 
to the PA of the outer disk. 

Such kinematic features, seen as extremes in HCO+, 
and probably also in CO(3-2), are expected in the frame¬ 
work of an additional radial velocity component directed 
towards the star. In order to impact the observed kine¬ 
matics, the additional radial velocity component must be 


comparable in m agnitude to the rotation velocity ( jCasas- 
sus et aI ^ |2Q13a ), so close to free-fall 


"Fig: m provides a summary of the data in each tran¬ 
sition. It is apparent that the coarser beam in band 7, 
compared to band 9, prevents a detailed study of the 
kinematics in the central 0.3". We also see in the decon¬ 
volved maps that attempts to trace the CO(3-2) kine¬ 
matics inside the cavity are thwarted by the absorption 
from the interve ning diffuse cloud in t he line of sight to 
HD 142527 (e.g. |Casassus et al.||2013b|). Once convolved 
with the clean beam the regions attected by the interven¬ 
ing CO cloud are less conspicuous, and may lea d to an 


accidental inter polation of the velocity field (as in Rosen- 


feld et aTlfeOldl). Due to this intervening absorption, the 

TO 


PA in CO(3-2) as a function of radius does not appear 
to rotate beyond 30-40° East of North, or ^60° relative 
to the outer disk. This shift is close to the value of 54° 
obtained when the radial velocity component scales as 
times the azimuthal velocity component (i.e. for a 
radial component close to free-fall). 

® with velocities corresponding to r’lsr/ km s“^ G [—2.8, 3.2] 
and G [4.8, 8.8], while the systemic velocity is at v\sj- = +3.8 ± 
0.1 km s“^ 


3. MODELING 
3.1. Parametric modeling 

In order to reproduce the intra-cavity kinematics, we 
developed on the accretion fl ows proposed to accoun t 
for the HCO+ observations (Casassus et al. 2013c), 
while also taking into ac count the warp imp lied by the 


scattered-light shadows (Marino et al. 


2015). We used 


a parametric model for HD 142527 that ht s the spectral 
energy distribution (see Marino et al.|2015, for a detailed 
description), with the addition of a radial velocity com¬ 
ponent corresponding to stellocentric accretion. Emer¬ 
gent specific intensity fields were predicted by transfer¬ 
ring the stellar radiation throu gh the model circumstel- 
lar material with the R ADMC 3package (version 0.39, 
Dullemond et al.1|2015 ). 


T’he warp is implemented by rotating the disk as 
a function of stellocentric radius. The orientation of 
the outer disk, with a fixed inclination i = 20° and 
PA=—20° beyond rorient_out = 23 AU, is linked linearly 
at Torient-in = 20 AU with a nou-coplauar inner disk, at a 
relative inclination of 70° and PA= —8°. The orientation 
of the inner disk is such that its inclination relative to 
the plane of the sky is i ^ 46.7°, so that the Eastern side 
is near, while East is far for the outer disl^p] The inner 
disk probably cont ains the orbit of th e lo w mass com 
panion reported by 


Biller et al. 


(2012) and Close et al. 


(2014). 

Radial cuts illustrating our choice of velocity and sur¬ 
face density profiles are shown in Eig. [^inside the cavity. 
We caution that this parametric model is meant to help 
the interpretation of the data, and is not derived from a 
systematic search in parameter space. It is not a dynam¬ 
ically consistent model. Along with the warp, we also 
add a radial velocity component that accounts for accre¬ 
tion. We connect zero accretion at ^140 AU to free-fall 
at 30 AU, with a cubic profile. The free fall halts at 
23.3 AU. The tilted inner disk is assumed to be in pure 
Keplerian rotation, so that inside 20 AU we set zero ac- 

http://www. ita.uni-heidelberg.de/ 

dullemond / software / radmc-3d / 

in other words the parallactic angle of the normal to the disks 
is rotated by 168° 
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Figure 2. Example CO(6-5) line profiles. We select 8 locations near the star to illustrate the line profiles obtained in the deconvolved 
datacube. The central plot shows CTmem, and is repeated from Fig. The spectra are extracted from single pixels in the deconvolved 
datacube. Solid black is used to plot specific intensity in 10^^ Jy sr“^, along the y—axis, as a function of LSR vel^ty in x—axis, in 
km s“^. Solid green plots the 1 a error on the specific intensity spectra, as estimated from Monte-Carlo runs (see Sec. [b|>. 


cretion (a build-up of accreted mass onto the inner disk 
could then trickle to feed pole-on stellar accretion). The 
total surface density profile in the gap follows from the 
assumption of steady state accretion: 


I](r) 


dM^ 

dt 


j{2TTrVr), 


as a function of spherical polar radius r, where w is the 
radial velocity component, 

= 0 if r < rgi, with rgi = 20 AU, or 
V'p — O.luff if r <C Tstop? with rg^op — 23.3 AU, or 
Vr = —0.9uff if r < rff, with r^ = 30 AU, or 


Vr = —0.9Uff X 


(frg2-r) 

(frg2-rgi) 


( 1 ) 


with rg 2 = 130 AU, and where / = 1.1 is a factor to 
avoid divergence at rg 2 , and 


v,= ,/2GM.xlj-j^ 


We use a stellar accre tion rate dM^/dt = 10“^ Mq (Gar- 
cia Lopez et al.|2006|). The stellar mass is set at 1.8 Mq; 
higher masses produce higher velocities inside the cavity 
than observed. 

The orientation of the disk is parametrised in spherical 
coordinates, where the vertical axis z' is the normal to 
the local disk plane, and where the origin of azimuth <p' is 
the local line of nodes (i.e. the disk PA as a function r). 
The azimuthal velocity is set to Keplerian rotation 


if r < rg 2 , 
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Figure 3. Line profiles in CO(6-5). Top: Observed profile, ex¬ 
tracted on the restored datacube in natural weights, using aper¬ 
ture photometry inside a 1.1'' radius (solid black), and also in 
the uvmem model with a 0.2" radius (in solid blue, magnified by 
a factor of 10 and offset for clarity). The horizontal dotted lines 
indicate the origins of flux density for each extraction. Bottom: 
Corresponding line profiles extracted on the radiative transfer pre¬ 
dictions described in Sec. The wiggles in the line wings have 
counterparts in the observations - this may be bette r appreciated 
by inspection of the channel maps in Appendix |A.1| 


V(f)' = vk^ except in regions where \vr\ > 0, where we 
assumed that the specific energy density was close to 0, 
so that V(p = (neglecting the internal energy). 

In the warp itself, so between rorientJn = 20 AU and 
I’orient-out = 23 AU, material continuously connects both 
orientations. We investigated a range of values for the 
vertical velocity component in the warp, and obtained 
a closer match to the observations with re lative ly fast 
velocities, = i;warp = r’Ksin(0') (see Sec. 3.2.2). 

The temperature profile we adopted is / = yO K if 
r < 30 AU, or 




-0.5 


K, if r > 30 AU. 


( 2 ) 


These temperatures are well below the dust temperatures 
in the closer-in regions. If we impose thermal equilib¬ 
rium between both solid and gas phases, the predicted 
flux densities soar an order of magnitude above the ob¬ 
servations. 

Close to the star CO is probably photo-dissociated. 
The line intensity of CO(6-5) appears to display a decre¬ 
ment inside 0.2", which we accounted for with a Gaus¬ 
sian taper of the gas abundance inside 30 AU, with ala 
width of 5 AU. 

3.2. Radiative transfer predietions 

3.2.1. Adopted abundance fields and depth of the cavity 

We adopted a nominal H 2 /^^CO ratio Xqq = 10“^, 
and the following molecular abundances relative to H 2 
(no selective photodissociation was required). For ^^CO 
we used the ISM isotopic abundance for ^^C/^^C of 76 
(Stahl et al. 2008), multiplied by Xco- Foi* C^^O, we pro- 
ceeded as for ^^CO, but with an ISM isotopologue a bnn- 
dance for C^^O/^^CO of 500 ( Wilson fc Rood]l9M ). For 
HCO+, we used a default value for the HCO+ abundance 
relative to H 2 of 10“^ (somewhat lower than observed in 
molecular clou d cores, and comparable to values inferre d 
in other disks, IWootten et al.||1979t jTeague et al.||2015|). 

After smoothing to the resolution of the deconvolved 
models. Fig. shows that the accreting warped disk 
model qualitatively reproduces the basic properties of 
the cavity seen in the available molecular tracers. The 
model is in agreement with the non-detect ion of ^^CO 


and C^^O inside t he gap. Both in band 6 ( [Perez et al. 
2015) and band 7 ( Fukagawa et al.|[2Q13 ), the ^'^CO and 
C^^O datacubes do not show any sign of either twisted 
kinematics, a central peak, or of the gap-crossing fila¬ 
ments. The CO isotopologue data display plain Keple- 
rian rotation, without signal inside ^90 AU. 

The decrements und er the crescent-shaped continuum, 
previously reported byjCasassus et al.|(|2013c|), are a nat¬ 
ural feature of the predictions. In con^tmuum-subtracted 
images, and with our gas temperature p rescription, the 
gas e ffectively acts as a cool foreground (Casassus et al. 


3.2.2. Central beam 

The warp inside the cavity of HD 142527, i.e. the 
continuous change in inclination at some specific radius, 
should significantly impact the observed kinematics. For 
instance, give n the orientation im plied by the scattered 
light shadows ( Marino et al.|2Q15 ), the blue-shifted flows 
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Figure 4. Moment maps in HCO+(4-3) and CO(3-2). For each specie, we plot both restored images (labeled rest .) and the correspo nding 
deconvolved models (labeled MEM). Columns correspond to spectral moments, from left to right we show: intensity (see Appendix |B.2| ), 
specific intensity peak, Gaussian velocity centroid, and Gaussian 1 cr velocity width. 


to the North should suddenly shift to red as the inclina¬ 
tion changes sign, at the radius where the far side be¬ 
comes the near side. Interestingly, as illustrated in finer 
detail in Fig. the observed centroid map indeed has 
structure very close to the central star, well inside the 
central beam. Within ^0.25", the ‘twist’ in PA, i.e. the 
rotation of the axis joining blue and red-shifted veloc¬ 
ities, appears to dwell at ^40°, and so offset by ^60° 
from the outer disk PA, as expected for a radial veloc¬ 
ity component that has reached close to free-fall values 
while also keeping in near-Keplerian rotation. Indeed, at 
a given stellocentric radius, in the presence of a radial 
velocity component that is times the azimuthal com¬ 
ponent, the peak line-of-sight velocity is offset from the 
disk PA by ^54.6°. However, inside 0.1", the PA joining 
blue and red shifted velocities appears to reach 110 ±10°. 

The extreme in PA rotation in the inner 0.1" cannot 
be accounted for by free-falling material that is also in 
close to Keplerian rotation. This can instead be under¬ 
stood as the halt of free-fall inside a warped structure, as 
the blue-shifted flows to the North-West (llh) suddenly 


doppler-fiip into red closer to the star, where the inclina¬ 
tion changes sign (far side turns near). After smoothing 
to the resolution of the deconvolved MEM datacubes, we 
obtain an approximate match to the velocity field for an 
abrupt warp connecting the inner disk at 20 AU with 
the orientation of the outer disk at 23 AU, as shown in 
Fig.[7fD. The warp orientation must align to that of the 
inner disk as far out as 20 AU, to account for the shift 
in declination at systemic velocities, which covers almost 
0.2" when crossing the stellar position. 

In a continuous warp material must flow from one plane 
to the other and be accelerated in the transition. So in¬ 
side a warp there should be a velocity component orthog¬ 
onal to the plane of the disk n^arp- This component is 
probably limited by the Keplerian velocity vk around 
the central star. For instance with = 0.1 Uff and 
'^0 = '^warp = in spherical coordinates, the warp’s 
radial width must be at least ^3AU to connect one ori¬ 
entation to the other. Interestingly, the inclusion of Uwarp 
approximates better the observed velocity centroid, with 
fairly point-symmetric high velocity ridges connecting 
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Figure 5. Radial cuts for the surface density profiles (red) and 
magnitude of the radial velocity component (blue). The CO abun¬ 
dance is modulated with a Gaussian taper inside 30 AU, to account 
for photo-dissociation. 



Figure 6. Intensity maps (moments 0) extracted from the 
RADMC3D predictions, after smoothing to the resolution of the 
deconvolved CO(6-5) datacubes. 

the cavity ‘twist’ with the inner disk, corresponding to 
the loci of bluest and redest velocities in Fig. [^b. These 
ridges are missing in the model without t;warp 


4. DISCUSSION 
4.1. Ro-vibrational CO 


The CRIRES observations by [Pontoppidan et ~aL 
(2011) assign bright rovib CO to an inner disk, with 
similar inclination as the outer disk but offset in PA, at 
61 ± 3° - so within uncertainties coincident with the di¬ 
rection of the high-velocity HCO^. The total line-widths 
in CO(i;=l-0) at 4.67/im, at zero-intensity, are close to 


^40kms“^ (Pontoppidan et al.|2011). With a spectroas 
trometric offset of only 0.2 AU, (JU[v=l-0) is either fairly 
uniform or else confined to the central few milli-arcsec. 

We check the astrometric signal of a modified version 
of our disk model and find it to be in agreement with the 
spectro-astrometric CRIRES observations. We take the 
predicted CO(6-5) emission and transform the intensity 
with an 1/r^ mask to mimic the CO ro-vibrational emis¬ 
sion, which traces higher temperatures and emits closer 
to the star compared to the CO rotational lines. We 
manually add continuum emission scaled to reproduce a 
line to continuum ratio of 1.1, as in the observed rovibra- 
tional CO, and convolve with a 0.192” EWHM Gaussian 
to mimic the PSE of the AO-assisted CRIRES observa¬ 
tions. After superimposing a 0.2” arcsecond wide slit 
positioned at 150 degrees, we recover a line profile and 
spectro-astrometric sign al similar in shape and m agni- 
tude to that reported by Pontoppidan et al. (2011) when 
using a ^ 2.5. 

4.2. What about an outflow? 

The orientation of the disk is well known. As explained 
in Secs. and Sec. |3.1[ the shadows seen in scattered 
light imply that the inclination of the outer disk is such 
that the far side lies to the East. The Eastern side being 
brigh ter in the mid-IR is a lso in support of this orienta¬ 
tion (Eujiwara et al. 2006). In addition several near-IR 
spirals~have been found t hat are all consistently trailing 


dEnkagawa et al.' |2006b| Casassus et al. 2012 |Canovas 
et al.'2013 ~ Avenhans et al. '2 014|), even in tn^ radi o 
extensions into the outer disk (jchristiaens et al. |2014). 
Given this orientation, we can rule out an outnowing 
disk wind: the velocity component orthogonal to the disk 
would broaden the lines and preserve reflection symme¬ 
try about the outer disk PA, while the radial component 
in the plane of the disk would twist the kinematics in the 
opposite direction than observed. 

4.3. Origin of the free-falling veloeities 

The radial velocity in a warped v iscous disk is pre¬ 
dicted to be greater than in a flat disk (Pringle 1992), but 
is not expected to exceed the sound speed Cg = [ti/r)vK- 
Supersonic disk inflow is in principle possible in re¬ 
gions of the disk that lose angular momentum rapidly 
through a magnetized wind, rather than via internal vis¬ 
cous transport, though whether this can occur without 
equally rapid mass loss (for which there is no evidence 
in HD 142527) is not clear. Alternatively, the free-fall 
velocity signal could be a consequence of strong dynam¬ 
ical perturbat ions. If the low-mass companion d etected 
in the system ( Biller et al.|2Q12 [Close et al.|2Q14[ ) is mis¬ 
aligned to the outer disk in the same way as The inner 
disk, which is known to be inclined by 70° (Marino et al. 
2015), then an observable dynamical impact is highly 
likely. The alignment of the companion in the plane of 
the inner disk seems a natural way to explain th e incli- 


nation of the inner disk itself. Miranda & Lai mm 


p 

]he 


have recently shown that large offsets between tne or- 
bital planes of circumprimary disks and the plane of the 
binary can be stable for several Myrs, longer than the 
expected lifetime of the disk. 

It has been shown recently that such strongly mis¬ 
aligned binary-disk systems can undergo a process 
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dubbed disk tearing dN ixon et al.| 2013 Nealon et al. 
2015 Dogan et al.|2015 ), where nodal precession torques 


Figure 7. Comparison of observed and model CO(6-5) kinematics in the central regions of HD 142527. The origin of coordinates is set to 
the stellar position. Velocity-integrated intensity in CO(6-5) is shown in grey scale. The colored contours for v\sj- have constant interval and 
are spread over [0.21, 7.87] km s“^ (as in Fig.jH. a): Observed moment maps extracted on the MEM datacubes. b): moments extracted on 
the radiative transfer prediction, after smootmng to the resolution of the MEM datacubes. c): same as b) but in native model resolutions, 
without smoothing. Regions without contours near the origin correspond to higher velocities, d): same as b) but with a slow velocity 
component perpendicular to the disk plane (uwarp in the text). 

The dynamical influence of the low-mass companion on 
its current orbit will not extend out to the 100 AU scale 
needed to explain the large extent of the cavity. It is 
possible, however, that its orbit may have changed over 
time. The highly inclined orbit of the companion relative 
to that of the massive outer disk c ould trigger Kozai os¬ 
cillations ( Teyssandier et aT]|2013[ ), with a period of the 
order of 10^ yr if the outer disk has a mass of 0.1 Mq. 
A Kozai cycle would result in (damped) oscillations be¬ 
tween the current inclined circular orbit, and a highly 
eccentric coplanar orbit which would perturb the disk 
out to greater radii. Numerical simulations are needed 
to determine if such a cycle is possible in HD 142527. 


induced by the binary produce a warp at the inner edge 
of the disk. If the disk is unable to communicate the 
warp efficiently, then the disk will effectively break at 
some radius inbreak- Annuli of gas then begin to tear 
off of the inner edge and can precess freely, undergo¬ 
ing self-interaction and angular momentum cancellation 
causing a sub sequent infall of g as driving very high ac¬ 
cretion rates (Nixon et al.||2012). 

To date this process has primarily been studied in the 
thin-disk regime appropriate for black hole accretion, 
where the dimensionless viscosity parame ter is larger 
than the d isk aspect ratio, i.e. a > HfR ([Nixon et al. 


2012 2013]). Nevertheless it has also been s hown to worE 
in th e thick-disk regime where a < H/R ( [Nealon et al. 


2015), which is more likely to hold true for circumstellar 
HTsks like HD 142527. 

It is possible to estimate a value for the radius inbreak 


at which d isk tearing will occur (see Eq. A3 of Nixon 
et al.[[2013 ) using the parameters of the disk and binary. 


With appropriate values for HD 142527, we expect that 
the disk could break at inbr eak ^ cl, where a is th e binary 
semimajor axis. However, [Nealon et al.[ ([2015[) found a 
consistent offset between the estimated minimum value 
and the true breaking radius in a series of SPH simu¬ 
lations, by a factor of 2-3. Estimates of a ^ 15 AU 
are therefore entirely consistent with a breaking radius 
of 30 AU, as demanded by our model. However, as this 
is likely very sensitive to the parameters of the system a 
more quantitative comparison of this scenario with the 
observed data is out of the scope of this paper, requiring 
the development of targeted numerical simulations. 

We also note that while this scenario seems to ade¬ 
quately explain both the inner warp and the free-falling 
velocities, it does not explain the presence of the ^ 
130 AU cavity in the system. This should not be sur¬ 
prising, given the large difference in scale between the 
warp and the cavity (and therefore different timescales 
associated with physical processes at these radii), but 
means that the model is not yet a complete description 
of the system. 


4.4. Long term Kozai oscilations and the large eavity 


4.5. Non-axial symmetry in the eavity. Signature of 
obseured planets? 

Observationally we cannot test directly for past copla¬ 
nar eccentric orbits, as in the Kozai oscilations, but we 
can instead look for other mechanisms that could account 
for the clearing, such as additi onal bodies at ^100 AU . 
Stringent limits are available (Casassus et al. 2013c), 
^4 Mjup bodies should have been detected unless they 
are obscured. 

An interesting feature of the Cycle 0 HCO+(4-3) data 
is the faint filament crossing the cavity to the East, at 
low velocites, with perhaps another filament to th e West 
(see Eig. [^for a summary, Casassus et al. 2013c). It is 
tempting to associate these HCO+ features with planet- 
induced accretion flows, i.e. gaseous streamers due to ob¬ 
scured protoplanets inside the cavity, that feed the faster 
flows at the disk breaking radius. As the critical density 
for the excitation of HCO+(4-3) is nH 2 ^ 10^ cm“^, and 
higher than the CO(3-2) critical density of ^ 10^ cm“^, 
in the streamer interpretation the HCO+ filaments and 
high-velocity flows would trace the densest regions of the 
general accretion flow across the gap, which is also re¬ 
flected in the ‘twisted’ CO kinematics. If the free-falling 
accretion is initially fed by planet accretion at ^ lOOAU, 
there should be corresponding azimuthal modulations of 
the H 2 density field. Perhaps this is reflected in the struc¬ 
ture of the faint, low-velocity HCO+(4-3) inside the cav¬ 
ity. 

An alternative to azimuthal modulations of the H 2 den¬ 
sity field in the cavity, is that only the abundance of 
HCO+ is modulated. The tilted inner disk is shadow- 
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Figure 8. Summary of Cycle 0 band 7 observations, from MEM 
maps, with continuum in red, HCO+(4-3) in green, and CO(3-2) 
in blue. Velocities have been restricted to highlight the fainter 
structures seen in HCO+, which are otherwise dwarfed by the fast 
HCO+ central emission. 

ing the UV radiation required to produce HCO+ in the 
cavity, resulting in HCO+ filaments with wide opening 
angles, as opposed to the thin protoplanetary stream¬ 
ers predicted in hydrodynamic simulations. A problem 
with this interpretation is that the outer ring is fairly 
round in HCO+, there are no obvious counterparts in 
HCO+ of the scattered-light shadows. Perhaps this re¬ 
flects a different production chain for HCO+ in the outer 
disk, where it could be driven by charge-exchange with 
cosmic-ray induced H^. 

5. CONCLUSION 

The new CO(6-5) data, along with the orientation of 
the disk inferred from the scattered light shadows, have 


allowed us to understand the intra-cavity kinematics in 
HD 142527. Stellocentric accretion, starting from the 
outer disk and reaching close to free-fall velocities, with 
a steady-state mass flow fixed at the observed stellar 
accretion, is consistent with the bulk properties of the 
available CO isotopologue data. The observed HCO+ 
flows are also consistent with this stellocentric accretion, 
but with an emissivity that is somehow modulated in az¬ 
imuth. Fine structure in CO(6-5) also suggests non-axial 
symmetry inside the cavity. 

While the data are consistent with a continous, yet 
abrupt and fast warp linking the two non-coplanar disks, 
further observations are required to understand the de¬ 
tailed structure of the intra-cavity kinematics and the 
inner warp. For instance, the large relative inclination 
observed between the inner and outer disks suggests that 
the fast accretion could be due to disk tearing. If the 
low-mass companion is contained in the inner disk and 
is breaking the disks, then its orbit would also be highly 
inclined with respect to the outer massive disk. The com¬ 
panion may then undergo Kozai oscillations, with high- 
eccentricity periods that may perhaps explain the large 
cavity. 
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tory, through program ALMA#2011.0.00465.S. ALMA is 
a partnership of ESO, NSF, NINS, NRG, NSC, ASIAA. 
The Joint ALMA Observatory is operated by ESO, 
AUI/NRAO and NAOJ. Einancial support was provided 
by Millennium Nucleus RC130007 (Chilean Ministry 
of Economy), and additionally by EONDECYT grants 
1130949, 1141175, 3140601, 3140634, 3140393. S.M. 
acknowledges CONICYT-PCHA / Magister Nacional / 
2014-22140628. PR and VM acknowledge support from 
CONICYT-ALMA grant alma-conicyt 31120003. AD ac¬ 
knowledges CONICYT-ALMA grant 31120007. PJA ac¬ 
knowledges support from NSE award AST 1313021. MM 
acknowledges CONICYT-Gemini grant 32130007. This 
work was partially supported by the Chilean supercom¬ 
puting infrastructure of the NLHPC (ECM-02). 


APPENDIX 


OBSERVATIONS 

ALMA CO(6-5) data 


Det ails on the instrumental setup are described in a companion article on the continuum emission (Casassus et al. 
2015). We used self-calibration to improve the dynamic range of the continuum images. Applying the same gain 
corrections to the line data also resulted in improved dynamic range. Before self calibration, the peak signal in the 
systemic velocity channel reached 0.42 Jy beam“^, in natural weights (beam of 0.30" x 0.25"), with a noise level of 
0.08 Jy beam“^ clearly dominated by systematics rather than thermal noise. After self calibration, the peak signal 
increased to 0.67 Jy beam“^, while the rms noise level decreased to 0.04 Jy beam“^. 

Continuum subtraction under CO(6-5) was performed with a first-order fit to the continuum in the visibility domain. 
The subset of channels neighboring the line was then split off into another datafile, and subsequently re-sampled in 
frequency into the LSRK frame. 

The CO(6-5) dataset is presented in channel maps in Eigs.j^and 
both the restored image and the underlying MEM model. 


10, with three-channel bins. We chose to present 
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Figure 9. Channel maps in the uvmem model datacube for CO (6-5). We binned the datacube into 3-channel averages. 

Alignment 

The comparison of the band 9 and band 7 continua reveals a positional offset of ^0.07 arcsec essentially in the 
East-West direction (precisely [-0.063,-0.016] in arcsec, along RA and DEC). As detailed in a companion article 
devoted to a multi-frequency analysis of the continuum, including a detailed presentation of the band 9 continuum, 
this positional offset is measured by cross-correlating the visibilities, but it is also noticeable by visual inspection of 
the continuum maps (for instance by blinking the images). With the 0.35 arcsec resolution provided by these Cycle 0 
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Figure 10. Channel maps in the restored CO(6-5) datacube. We binned the datacube into 3-channel averages. 


observations in band 7, the center of symmetry of the CO(3-2) velocity centroids coincide fairly well with the position 
of the central star at the epoch of observations. Hence for the multi-frequency continuum studies we chose to take the 
band 7 astrometry as reference. 

However, at 0.02 arcsec resolutions the center of symmetry of the velocity centroid appears to be offset from the 
nominal stellar position as given by the default band 7 world-coordinate-system (WCS). This represents an opportunity 
to improve the astrometry by tying the stellar position with the center of symmetry. We aligned the velocity centroids 
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by applying an offset of [—0.109, —0.031] arcsec to the default band 9 WCS. As illustrated in Figs. Hand 11, the 
center of symmetry of the velocity centroid coincides with the origin of coordinates, located at the st^ar position. 
The peak velocity dispersion is also close to the origin, but in this case the exact location of the peak depends on the 
moment extraction algorithm, as the signal is spread over a large range in velocity and includes low signal-to-noise 
channels. 


IMAGE SYNTHESIS 

Angular resolution and statistics of the deconvolved MEM models 


Details o n our image synthesis algorithm uvmem are given in a companion paper on the continuum emission (Casassus 
et al.|2Q15 ). Here we exploit the finer angular resolutions provided by the MEM models, which requires an assessment 
of their statistical properties. 

We start by scaling the observed visibility weights so that they correspond to the observed scatter. In order to 
preserve the relative values of the calibrated weights, a constant scale factor is applied to all visibilities. This factor is 
the median of the ratio between the original weights and the dispersion values. 

Estimates of the angular resolution of the deconvolved images were obtained as follows. We selected a reference 
channel, at the systemic velocity, and split off the corresponding visibilities into a reference visibility dataset. We 
then estimated the PSE by simulating interferometer observations on a point source, with an identical v coverage 
as the reference visibilities, and with the addition of Gaussian noise, with a dispersion given by the visibility weights. 
In practice the point source is a spike, defined on the same canvas as used for the deconvolved images, and whose 
fiux density corresponds to that of the whole observed image. A deconvolved model of these simulated point-source 
visibilities is then obtained using the same algorithm as for the observations (i.e. uvmem). This deconvolved model 
represents the PSE of the deconvolution algorithm, with a central source and secondary maxima fainter by a factor 
of at least 20. An elliptical Gaussian fit to the central peak defines the finest resolution achievable in the deconvolved 
images, of 0.062" x 0.051". 

Since the maximum-entropy regularization is non-linear, the uvmem PSE depends on the model. We approximate 
the case of extended emission by spreading the total fiux in a collection of spikes. The ‘clean’ beam in this case varies 
by up to 50%. So, for the purpose of comparing synthetic radiative transfer predictions with the uvmem model images, 
we chose a representative resolution of 30% times coarser than the finest PSE. 

Monte-Carlo simulations provide estimates of the uncertainties associated with these MEM model datacubes. We 
added Gaussian noise to the observed visibilities, as given by the visibility weights, and obtained a deconvolved model 
image with uvmem, with the same settings as for the observations. We then took statistics on 120 realizations, checking 
that the average datacube satisfactorily converges on the original model datacube without the addition of noise. 


Moment maps 

In the case of noisy data, the standard approach to extract moment maps by collapsing the velocity axis injects 
unnecessary noise when the line-of-sight line profiles are much narrower than the spatially integrated line. Rather than 
collapse the datacubes with standard moments, we extracted moments by restricting the velocity range above a given 
threshold (which is set to 0 for the restored maps, with thermal noise). Eor each line of sight, we identify the position 
of the peak specific intensity. We then search for the first nulls, both towards the red and blue, and extract a spectrum 
leaving room for a few more channels (two in this case). Last, we take standard moments over this reduced velocity 
range, and also fit a Gaussian. The resulting moment maps are embellished by median-smoothing with a square 5x5 
kernel. 

In the case of the MEM deconvolved models, the strict positivity requirement forces an intensity floor /min, which 
we take at 1/100 the estimate noise level (in natural weights). We subtracted a-posteriori /min from the deconvolved 
maps. However, the positivity requirement biases the noise, so that velocity dispersions are artificially broad in region 
of low signal (such as the outer ring in HD 142527 in CO(6-5)). To account for this we applied the velocity restrictions 
described above, with a threshold defined by the rms dispersion of the noise in velocity channels devoid of signal. 

Comparison of different image synthesis schemes 

The moments maps extracted on the clean, and uvmem datacubes are compared in Eig. EH The details of the 
centroid map do not depend on the reconstruction. However, the intensity and velocity dispersion maps inside the 
central beam (so within 0.2" of the star) result from the superposition of several velocity components. 


IMPACT OE A EINITE BEAM IN ESTIMATES OF INTRINSIC LINE-WIDTHS. 

An alternative cross check on the quality of the deconvolved velocity dispersion as an estimate of the intrinsic line- 
width, at least away from the central beam, is to take statistics on the velocity centroid i;mem extracted on the MEM 
datacubes. We degrade i;mem with the same beam as used in the Briggs restoration, — ('^mem) = '^mem G H, and 

calculate the velocity dispersion in the restored maps due to the bulk flow with dbuik = \/('^mem) “ ('^mem)^? where 
('^mem) — '^mem G B. Wc then subtract (Jbuik in quadrature to estimate the intrinsic line-width: dint. = \/crrest. ~ crbuik, 
where drest. is the 1 a velocity dispersion extracted from the restored datacubes. The corresponding dispersion at P 
is crint.(^) = 1-81 km s~^. A few pixels have slightly lower drest. than dbuik, which may stem from the finite resolution 
of i;mem, or to nou-Gaussiau line profiles. As these pixels represent a small fraction, we set them to zero. Eig. 
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Figure 11. Moment maps in CO(6-5) emission, x— and y—axis show angular offset along RA and DEC from the stellar position. The 
upper row shows maps extracted on the clean datacubes, produced in Briggs weights. The middle and bottom row respectively show 
the maps extracted on the restored uvmem+CASA datacubes, in Briggs weights, and the deconvolved MEM models. Columns correspond 
to intensity (left), velocity centroid (i nidd leb and 1 cr velocity dispersion (right). The moments from the deconvolved MEM maps are 
extracted without the thresholds (Sec. |B.^ used in Fig. ^ 

compares the deconvolved velocity dispersion, also shown in Fig. with the dispersion of the MEM datacube after 
correction for bulk motion. 


The (Jrest. iiiap in Fig. shows interesting structure on scales much smaller than the clean beam. Two pe aks in 
velocity dispersion are very r eminiscent of the orientation and separation between HD 142527 and HD 142527B (Biller 


velocity dispersion are very r eminiscent or tue orientation and separation between riu l4ZbZi and riu (Diller 

et al.|2Q12 (Close et al.|2Q14| ). The binary separation is 84 mas along PA=127°, and the mass ratio is < 0.1, where the 


inequality is due to the correction for accretion luminosity on HD 142527B. Although it is tempting to associate the 
secondary peak in CO(6-5) velocity dispersion with HD 142527B, we see that it does not persist in the deconvolved 
maps, while the main peak also shifts slightly in position. These peaks must reflect features of the underlying velocity 
field, so that the corresponding pixels mix blue and red shifted emission in a coarse beam. 

Facilities: ALMA. 
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